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Abstract. We consider the compressible Navier-Stokes equations for isentropic dynamics 
with real viscosity on a bounded interval. In the case of boundary data defining an admis- 
sible shock wave for the corresponding unviscous hyperboUc system, we determine a scalar 
differential equation describing the motion of the internal transition layer. In particular, for 
viscosity e small, the velocity of the motion is exponentially small. The approach is based on 
the construction of a 1-parameter manifold of approximate solutions and on an appropriate 
projection of the evolution of the complete Navier-Stokes system towards such manifold. 
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1. Introduction 

This article is devoted to the description of slow-motion for the hyperbolic-parabolic Navier- 
Stokes system for compressible isentropic fluid with real viscosity, that is, in term of the variables 
density/velocity (p, w) 

dp dipw) „ dipw) 5 f 9 5wl „ 

The functions definining the pressure P and the viscosity v are required to satisfy the standard as- 
sumptions 

P'ip) > 0, P"(p) > 0, v(p) > 

for any p under consideration. The relevant cases P(p) = C p" with a > 1 and y(p) - C pP with 
a,p > I and C > fit into the general framework In particular, shallow water Saint- Venant system 
with viscosity corresponds to the case a - 2, ^ = \ (see [8] for the derivation of the model and ||5l 
for a recent review on shallow water equations). 

Given € > 0, the space variable x belongs to the bounded interval {-€, () and system ([T) is comple- 
mented with boundary conditions 

p{-€) = p„ w{±i) - w± > 0. 
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Let us set 

F(u,v) = — +P(u). 
u 

Considering the variables density /momentum {u, v) - {p,p w), system ([T) becomes 
with the boundary conditions 

(3) w+ u{+i) - v{+i) = 0, v{-{) = p^w^ 

In the limiting case ( +oo, system Q is known to support traveling wave solutions, i.e. solution 
with the form {u,v) - {U{x - ct), V{x - ct)) satisfying the asymptotic conditions 

([/, V){-oo) = v_), {U, V){+oo) ^ {u^, v+) 

for appropriate choices of (m±, v+) and c given by the usual Rankine-Hugoniot relation. Thanks to the 
galileian invariance of Q, we may assume, without loss of generality, the speed c to be zero. In such 
a case, the component V turns to be constant, so that v+ and v- are forced to be equal to a common 
value, denoted here by v*. Hence, the equivalent form of Q for the special solution (U, V) is 

U{±{) = u±:=—, V{-€) = Vh. 
w± 

The wave solution {U, V) belongs to a 1 -parameter family of stationary solutions to Q generated by 
the space translation group. The stability analysis of such family has been explored for a long time 
and a number of orbital stability results for different regimes and structure functions has been proved 
(see lIlOllIlITOlEIl). 

The analysis of the dynamics in bounded domains has been also investigated (among others, we 
quote ll4l[T6l[l3)> but always with a Mmited attention to the dynamical behavior close to equilibrium 
configurations. In term of special solution, in bounded intervals {-i, £) with fixed £ > 0, the situation 
is different with respect to case of the whole real line. Formally setting e = in (|2l), there still 
exists a one-parameter family of steady states given by a sharp transition at some point ^ € (-{, {). 
Differently, for e > 0, there exists a single steady state satisfying the boundary conditions Q with 
w±,v- (or, equivalently, m+,v*) chosen so that there exists a stationary wave in the whole real line 
with same asymptotic states. 

As a consequence, in analogy to the case of scalar viscous conservation laws, it is expected that, 
in the regime of e small, the solution determined by an initial datum consisting of a single transition 
from (m_, V*) to (m+, v*) converges in a short timescale to a specific profile with transition located at 
some point ^ - ^(t) and, then, on a much longer timescale, moves to the location of the single steady 
state. Such motion is an effect of the boundary data and it is expected to be very slow. Our aim is to 
derive formally a differential equation for the location ^ and to show that the motion of ^ is indeed 
exponentially slow for e ~ 0. The first tool is the definition of an approximate invariant manifold 
{W{-,^) : ^ e (-^,0} whose elements are approximate steady state of © and resemble, in a sense, 
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transitions from v*) to (u+, v*) located at ^. Among the many different and significant choices for 
constructing such manifold, our preference goes to the one we explored in |fT9l for the case of scalar 
conservation laws and that consists in matching exact steady states in (-{, ^) and in {) at x = ^ 
by imposing some appropriate conditions. Assuming that the spectrum of the linear operator 
obtained by linearing Q at W{-,^), has a first eigenvalue that is real and simple, applying a projection 
method in the spirit of ||23]| . we determine the equation (see Section|3]l 

dt ^) M+ - U- 

where {(p{-,^),tfr{-,^)) is the first eigenfunction of the adjoint operator and the functions k± are 
implicitly defined by 

— ds = ± 

Ju, s'^{k+ - F{s,v^,)) E 

where is the value determined by the condition P'{u^,)ul = v1. The difference function ^ i-^ 
/c+(^) - ic-{^) is continuous, monotone increasing, diverges at ±00 as ^ ^ ±i; hence, the above 
differential equation possesses a single equilibrium point corresponding to the unique stationary 
solution of the problem. 

For small s, the elements of the approximate manifolds tend to a piecewise constant configuration, 
with a single jump located at ^. Thus, it is possible to determine the leading term in the expression for 
the eigenfuntion (0, i/') of the adjoint operator and to obtain a new version of the motion equation 

d^ _ {^+i)AmK^{^)-K-{0 

dt duF{u-,v^:) u+ - u- 
where is the first eigenvalue of the operator Xf. Let us stress that, at such step, the choice 
of the boundary conditions is particularly relevant, since it determines the specific structure of the 
eigenfunction {cj), ifj). 

As it should be, the behavior of A\ for small s plays a crucial role (details are given in Section 
O. Extrapolating from llT4l . that concerns with general systems of conservation laws with a (non- 
physical) second-order parabolic term, and by the quoted results on asymptotic stability of steady 
waves on the whole real line, we expect that A\ is (exponentially) small and negative. We are not 
aware of any availabe rigorous result in this sense at the present time. 

Assuming that such eigenvalue stayes bounded in the limit e ^ 0^, the speed of motion along 
the approximate manifold, is determined by the functions k+ , whose leading terms can be determined 
starting from their definitions, obtaining the final form for the motion's equation 

dt yOur - U- \ Ou(j £ 1 U+ - \Ou(j S J) 

where duG{u, v) = vv{u)/u^ and the upper scripts - indicate that the function is calculated at (m+, v*). 
In particular, for small e, the motion is exponentially slow. The stability of the equilibrium point is 
still encoded in the sign of the first eigenvalue Ai that the present analysis is not able to reveal. 



4 



This paper follows the research line on metastable behaviors for conservation laws widely explored 
in the last decades. The first contribution has been the pioneering article [131 concerning the analysis 
of the scalar Burgers equation, that has been also the subject of |23| (based on the use of projection 
method and WKB expansions) and [15 J (standing on an adapted version of the method of matched 
asymptotics expansion). A rigorous analysis has been performed in ||6] |7]), where the one-parameter 
family of reference functions is chosen as a family of traveling wave solutions to the viscous equation 
satisfying the boundary conditions and with non-zero velocity. Slow motion for the viscous Burg- 
ers equation in unbounded domains has been also considered in literature: the case of the half-line 
(0, -i-cxj) has been treated in Il24l [T8l [221 : while the case of whole real line has been examined in 
|[T2l [TTl 121 (with emphasis on the generation of A/^-wave like structures and their evolution towards 
nonlinear diffusion waves). 

Despite of the wide number of contributions to the stability of traveling profiles in the whole real 
line, results relative to slow motion and metastable behavior in the case of systems of conservation 
laws appear to be rare. We are only aware of |l9l (that uses asymptotic expansions to deal with systems 
of conservation laws, with model examples being the Navier-Stokes equations of compressible vis- 
cous heat conductive fluid and the Keyfitz-Kranzer system, arising in elasticity), |[T4l (that deals with 
the problem of proving convergence to a stationary solution for a system of conservation laws with 
viscosity, with an approach based on the analysis of the linearized operator at the steady state) and 
El (that addresses to the Saint- Venant equations for shallow water and, precisely, the phenomenon of 
formation of roll-waves, by means of a combination of analytical techniques and numerical results). 

In this respect, we consider our contribution, even if mainly based on formal arguments, original 
and hopefully stimulating for people working in the area of dynamical properties of solutions to 
systems of conservation laws. Specifically, it seems that the equation we propose for the motion of 
the transition layer ^ (that is the slow dynamics along the approximate equilibrium manifold) is the 
first attempt in this direction for isentropic Navier-Stokes equations for compressible fluids. 

The article is divided into three more Sections. In Section 2, we present the general procedure 
to derive formal equation for the motion along an approximate equilibirum manifold in the case of 
general hyperbolic-parabolic systems. For pedagogical reasons, we also show how the approach 
simplifies in the case of scalar Burgers equation. Section 3 is the heart of the paper. It deals with 
isentropic Navier-Stokes equation for compressible fluids, written in the form ((2]i. After recasting 
the admissibility condition for entropic jumps in the unviscous case, we build up the approximate 
equilibrium manifold, working with time-independent solutions and matching them at x = ^ hy 
means of appropriate transmission condition. Using such special approximate solutions, we are able 
to determine an equation describing the evolution along the manifold. In order to determine an explicit 
expression for the ratio i/f(^,^)/^(^,^) appearing in such equation for the motion along the manifold, 
we analyze the eigenvalue problem for the adjoint operator and we deduce an approximated version 
of its solutions for the regime s small, by approximating the element of the manifold to be piecewise 
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constant functions with a single jump located at ^. Finally, in order to get an ultimate version of the 
motion's equation, we determine the leading term for the functions k+ in the limit e ^ 0"*^, Section 4 
contains the conclusions and it is mainly dedicated to propose a number of eventual research direction 
motivated by the present work. 

2. Scalar reduced dynamics for hyperbolic-parabolic systems 

As a first step, we present the strategy to obtain approximate equation for the motion along an 
approximate manifold of solutions to a general hyperbolic-parabolic system having the form 

(4) _,_|/,„,_,Kw,-| = 0. 

where w - w{x, t) e R", x e /, ? > 0. Our target is to apply such approach to the case of isentropic 
Navier-Stokes equation. Thus, for the sake of simplicity, we do not state precise assumptions on 
the structure functions / and b and we proceed in a purely formal way. In any case, we expect the 
procedure to be meaningful for the usual class of hyperbolic-parabolic systems considered in the 
recent literature (see ||25| and descendants). 

Given T > 0, we consider the initial-boundary value problem for @ determined by the conditions 

(5) w|^^(, = Wo 

complemented with appropriate boundary conditions. 

Given an open interval J and a one-parameter family of functions 

satisfying the boundary conditions, let ^ i-^ R^{- ; ^) be the distribution-valued map defined by 

(6) {R{-,^),if) := - ^j^Eb{W)—-f{W)]-—dx 

for any continuously differentiable function O : / — > R" with 0(±^) = 0. In what follows, we call 
R{-;^) the residual of W{-;^) with respect to equation The family is considered as an 

approximate invariant manifold for in the sense that the residuals /?(•; ^) vanishes as e ^ 0"*^ (in a 
sense to be made precise) for any ^ e 7. 

Next, we look for solutions to dUl-© in the form 

wi;t) = Wi-,m) + zi;t) 

with unknown ^ = ^{t) and z - z.{-, t) to be determined. Substituting into (H)) and disregarding the 
nonlinear terms in v, we obtain an approximated equation for the perturbation z 

(7) _.x,z + /?(-,^)-— (s^)- 

where /?(•; ^) is the distribution defined in Q and Ji^ is the linearized operator at W(-; ^), i.e. 

d { dz 
L^z := — shm^ + dhm z — - df{W)z 
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where we use the notation 

WVwV- ■= V 

dwk 



{db{W)zw)i ■= J] -^(^^^kWj with b = (bij). 

j,k 



Next, let us assume that, for any ^ e /, the operator has a first eigenvalue A\(^) that is real and 
simple. Let be the adjoint operator of 



where ' denote the transpose and 



Given ^ > 0, set / :^ {-£, i) and 

{u,v) := J u{x)-v{x)dx v € [L^(/)]", 

where • denotes the usual scalar product in R". Denoting by wi = a>i(-;^) an eigenfunction of the 
adjoint operator relative to the first eigenvalue and setting 

zi -zi(^;0:-(wi(-;az(-,0>, 

we determine a differential equation for the function t i-^ ^{t) by imposing that the component z\ is 
identically zero, that is 

^{ujy{-,m,z{;t)) = and <wi(-;^o),z(-,0))> - 0. 
at 

Using equation (O, we infer 

Since {u\,£.^z) - A\{a}\,z), we obtain a scalar differential equation for the variable ^, the latter 
equation can be rewritten as 



{, dW^ doji 1 d^ , „^ 



to be considered together with the condition on the initial datum 

(wi(-;^o),z(-,0)> -0. 

Neclecting the second term in the coefficient of the derivative of ^, we end up with an (approximated) 
equation for the motion along the manifold 

dVf d£ 

(8) <wi(-;a ^(•;^)>4 = {m{-,0,R{-,0)- 

ot, dt 

Our next effort is to determine a practical version of the above equation in the limiting regime e ^ 0^ 
and for the specific case of hyperbolic-parabolic systems ©. 



First of all, given ^ € I = (-£, i), we require the function W{-;^) to converge in the limit e — > 0"*" 
to the step function jumping at x = ^ from w_ to w+ 

(9) hm W{x, ^) ^ w-X(_ceM) + w+Xi(+cM) 

e— >0+ 

If the limit is in the sense of distributions, approximating the increment ratio for /j > by 

W{x-^ + h)-W{x,^) 1 ^ 
1 -/^M^(«..,(-) 

we infer the asymptotic representation 

dW 

—(•;^) = -[w]<5f (0 + 0(1) ase^O+, 

where [w] := w+ - w_ and 5^ is the Dirac distribution concentrated at ^. Thus, for small s, there hold 

dW 

Finally, assuming W to be chosen so that /? is a Dirac distribution concentrated at ^, 

(10) R{x,0-m6^{x) 

for some function r = r(^), we deduce our final expression for the reduced dynamics (|8) along the 
manifold {W{-;0] 

(11) -| = g(^):=- 

In order to make equation ([TTI ). it is necessary to determine the specific expression of the function r, 
describing the residual of W{-;^), and the ratio between the components of the first eigenfunction a>\ 
of the adjoint operator in the direction of r and in the direction of the jump [w]. 

Viscous Burgers equation. Let us consider the case of the scalar Burgers equation with viscosity 

(12) ^ + ^{-w^ -e — ^ = 0. 

at ox yi OX) 

with boundary conditions 

w{-i, t) = w, w(+£, t) = -w 
for some given w > 0. In this case [w] = -2w and equation ([TT]) reduces to 

(13) ^ = 

dt 2w 

where r(^) is the residual of W{-;^). 

In order to satisfy the requirement (fTOl i. we consider a specific approximate invariant manifold: 
for ^ e /, we build W{x, ^) by matching two steady states of the equation with appropriate boundary 
conditions. Namely, we set 

lW-{x,^) -i<x<^<i 
\W+{x,^) -t<^<x<t. 



where W± are steady state of (fT2]) in (-£, ^) and {), such that 

W.{-t, ^) - w, W-(^; = ^) = 0, ^) = -w. 

Functions W± can be expressed by means of an implicit formula. For w > 0, let us set 

X :- {(w,/<;) : w € (-w, vv), /c > w^/2j . 
Then, defining the function F = T{w, k) with (w, /c) e X by 



r(,„):= JI,a„h-(-^] 

Jo K- V \ V2^/ 



functions VK± are implictly given by 

(14) EnW^ix,aK±)=^-x. 
where the values k± = k+{^) are uniquely determined by the conditions 

(15) EY{±W,K^) = ^±i. 

Since W- and W+ are steady states of (fT2l) in {-i, ^) and (^, ^), respectively, we deduce, integrating 
by parts, that the residual R is 



dx 2 dx 



-£^(^)cD(^) + s^i^mo = £ 



dx dx 
Differentiating ([T4] | with respect to a;, we get 

(16) 



dW 
dx 



O(^) 



dW+ 1 7 



thus, in the notation of (fTOl) . there holds 

[5W 



r(^) - e 



5x 



dW- 

ox ox 



giving an "almost explicit" expression for ([T3T i. 

Since we are considering the regime e ^ 0^, it is possible to approximate the formulas defining 
K± and obtain a simpler o.d.e. describing the slow dynamics. Handling the explicit expression for the 
function F, conditions (flSl) become 



2k+ tanh 



- w. 



In particular, in the limit e — > 0"*", we get 



lim /f+(^) = -w . 

s—>0+ Z 



Thus, as e ^ 0^, there approximately holds 



2tanh{w(^±^)/2e 



_ « iw2|l+4exp(-^(^±^) 
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Collecting, we end up with the equation 



^ «w|exp|-% + ^)J-exp|-^(^-a 
corresponding to the formula determined in [,23 ( . 

3. Compressible isentropic Navier-Stokes equations 

Given smooth functions P - P{u) and v - v{u), let us consider the hyperbolic-parabolic system for 
compressible isentropic fluids 

. du dv ^ dv d iv^ , , , d lv\\ ^ 

(17) — + — = 0, + \- + P{u)-Ev{u)—(-]\ = 0. 
ot ox ot ox [ u ax\ulj 

where the pressure P is such that P'{u) > and P"{u) > and the viscosity v is such that v(m) > for 
any u under consideration. 

System (fTTl) is considered for x € {-£, €) together with the boundary conditions 

(18) v±u{±€)-u±v{±€) = 0, v{-€) = v- 
for some u±, v± with the value v_ being strictly positive. 

Admissible jumps in the vanishing viscosity limit. As a first step, let us consider the limiting 
unviscous regime £ ^ 0^. Putting formally £ = in ([TT] ). we obtain the hyperbolic system for 
unviscous isentropic fluids 

(19) du ^ dv ^ dv ^ d iv^ _^ p^^ ^\ q 

dt dx ' dt dx\u J 

Physical solutions are weak solutions to ( [T9l l satisfying appropriate jump conditions determined by 
the couple entropy/entropy flux 

2 3 

£(m, v) :^ ^ + n(M), Q{u, v) := At + I\'{u)v. 
2u 2u^ 

where IT is such that Y\"{u) = P'{u)/u. In particular, given m+ > 0, v+ e R and c € R, the function 

(20) (m, v)(x, t) := v_)^,_„„/x) + v+)^(„^„/x) 

(where is the characteristic function of the set /) is an admissible solution to (fT9l ) if it is a weak 
solution and it satisfies the inequality 

(21) dtE + d^Q < 

in the sense of distributions. Thanks to the galileian invariance, we may assume without loss of 
generality, the speed c to be zero. Thus, the requirement of being a weak solution translates in the 
classical Rankine-Hugoniot conditions, that read as 

,,2 



(22) [v]o - 0, 



V 

— + P{u) 
u 



= 
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where [g]^ := g+ - g- denotes the jump of the function g at x - ^. Similarly, the entropy condition 
(I2T]) reads as 

1.3 , 1 

— + n'(M)v < 0. 

-M Jo 

Conditions (|22]) -(|23]) select the possible couples (w-,w+) such that jump solution (l20l ) defines an 
admissible solution to ( fT91 ). The admissible couples can be explicitly characterized. 



(23) [Qh 



Lemma 3.1. Given P e C^([0, +00)) such that P'{u), P"{u) > Ofor any u>0, let us set 

(24) F{u, v) := — + P{u) m > 0, v € R. 

u 

Then, for any vi > 0, there exists /(vi) > such that for any V2 > f{v\) equation F{u,v\) - V2 
possesses exactly two solutions u± :- m±(vi,V2). Moreover, the function f is strictly increasing and 



convex and such that /(O) = P(0), /'(O) = 2 
Proof. There hold 

dF v^ , d^F 2v^ 

—(u, v) = --+ P'iu), —iu, V) = — + P"iu) 

OU OW^ U^ 

Hence, for any v > 0, the function F{-,v) has a single absolute minimum point m* = Ut(v) uniquely 
defined by the implicit relation 

P'iu,)ul = v^. 

Then, the function / = f(v) is defined by 

/(v) := minF(M, v) = — + P(m*) - P {u^)u^ + /"(««)■ 

«>0 Uif 

Moreover, differentiating we deduce 



df , ,, , , dUf 2v 

-i- = {P"{u,)u, + 2P'{u,)] — 1 = — > 
dv dv Ut 



and 

d^f 2P" 



dv^ P"u^ + 2P' 
showing the stated properties of the function (p. 



> 0, 



Proposition 3.2. Let P e C^([0, +00)) be such that P'{u),P"{u) > Q for any u > 0. For any couple 
(vi, V2) with vi > and V2 > /(vi) there exists unique {u+, v±) such that 

^\ 

V+=V\, F{u±,Vi) = ^ P{U±) = V2, U- < u+ 

u± 

such that the function (1201 ) is a weak solution to ( 1191 ) satisfying condition (I23I) . 

Proof. Couples {u±,v±) connected by a single entropic stationary jump are such that v_ and v+ are 
equal with common value denoted by v* and 

F(m_, V*) - F{u+,Vt). 
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The properties of the function F shows that, given there is a single couple of values for which the 
above relation is satisfied. It only remains to analyze condition (|23] ). 
For vi > 0, the entropy condition becomes 

2 2 

(25) A(m+, vi) := -^^ + 2n'(M+) <^■ + 2IV\u-) = A(m_, vi). 



Since there holds 



5A 2v\ „ 2dF 

— (m, vi) = + 2n (m) = -— (m, vi), 

OU u ou 



the function A has the same monotonicity of F with growth rate that decreases when u increases with 
respect to the one of F. Thus, condition ( |25] ) is satisfied if and only if m+ > m^. □ 



As an example, in the power-law case P{u) = ku""^^ /{a + 1) with a > 0, there holds 

1 2 a + 2 _i_ 2(g+i) 

a + 1 

Note that the function / is not differentiable two times at m = 0. 

Approximate invariant manifold. Next, we build a one-parameter family of functions W = W{-;^) 
for ^ € {-£, £) forming an approximate invariant manifold for dTTi and converging as £ ^ to 

(26) W'hyp(-^' ^) = (M-^v^)X(_gj^{x) + (m+, Vh.)^(^_,,)(x) 

Recalling the general procedure presented in Section |2] we want to choose W so that the residual R{^) 
is a delta distribution concentrated at ^. Thus, given ^ € /, we opt for defining W by matching at ^ € / 
at the state (m*, v*) the two stationary solutions of ([TT] ) in {-£,^) and {^,-£). Precisely, we denote by 
W- - {U-, V-) and W+ = {U+, V+) the solutions to 

dv d (v^ d /v\] 

(27) — -0, \-+P(u)-sy{u)—{-)\ = 0, 

dx dx yu dx\ul] 

in {-i, ^) and in (^, {) respectively, satisfying the boundary conditions 

W.{-€;^) = (M_,v,), W-(^;^) - (m.,v,), 

W^i^-e = (M*,v,), W+(^;^) - (M+,v,), 

where is the absolute minimum point of the function F{-,Vf), i.e. is uniquely determined by the 
requirement P'{ut)ul = vi (see Lemma [3TT]) . 
Solutions to dTTl) solve 



dG du 

(28) v = vi, e— (m,vi)— = /^-f(M,vi) 



where v\,k are integration constants, the function F is defined in (I24b and the function G is given by 

f" v{s) 
G{u,v) :- V — — ds, 
Ju 

for some fixed m > 0. From now on, let us consider vi > 0. 
For F defined in (l24l ) and given v* , let us set 

S:={(m, /c) : u€{u-,u+), a;>F(m, v*)} 
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and define the function 



T{u,K) := — 
Ju. K- 



G{s, v*) 



ds 



for (u, k) e E. 



Tlien, solutions to (ITT] ) satisfying the condition m(^) = are implicitly defined by 

V(x) = Vt, £r{u{x),K) = X - ^ 

Denoting by u±{ic) the two solution of F{u, v*) = k (see Lemma [37T]) . the function F is such that 
Y{u, +00) = 0, Y{u-{k),k) = -00, T(u+{k),k) = +00 

T{u, •) is increasing if m < r(u, •) is decreasing if m > m*. 

As a consequence, for any ^ € (-^, there exist (unique) k± = k±{^) e (f (m+, v*), +00) such that 

eY{u-,k-) = -£ - ^ and er(M+, a:+) = € - ^. 

Correspondingly, we set 

({U-(x,^),v,) -€<x<^<e 
Wix n = { 

\{u^{x,^),v,) -e<^<x<e, 

where functions ?7+ are implicitly given by 

(29) sY{U^{x,^),K^) = x-^. 

Next, let us calculate the residual of W: for any test function <^ = {(p\,(p2), there holds 

r f dG dU 1 d(fi2 

(/?(•; a 3>> - {FiU,v,)-s—iU,v,)—]^dx 



dG dU^ 

F([/_,V,)) -£ — ([/_, V,) — 
OM OX 



dG dU+ 

OM OX 



-H-(^^'V,)— 
ou ox 



Differentiating (1291 ) with respect to x we infer 



£^(t/±(x,f),v.)^(^,^) = ^± -F([/±(x,^),v.) 
ou ox 



so that, with r = r(^) as in (fTOl ). we obtain 

r(^) = (0,^+(^)-^_(^)) 

The functions ^ i-> /c+(^) are the inverse of the relations 

^ - -£ - eT{u-,k-) and ^ - £ - sT{u+,k+), 

thus, as a consequence of the properties of function F, the difference function ^ 1-^ /<'+(^) - K-i^) is 
monotone increasing and such that 



lirn^ K+(^) - K-{^) = -00 and 



hm K+{^) - K-{^) = +00. 



Therefore, there exists unique € {-£,£) such that (/c+ - K-){^t) = and such a value is such that 
(f/^(-; ^H.), Vh.) is the unique steady state of the problem. 
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Finally, denoting by coi the first eigenfunction of the adjoint operator and setting coi(x,^) 
((p(x, i^(x, ^)), equation ([TT]) becomes 

_ a ipi^, ^)) • (0, K^io - K.{0) 

dt 

that gives 



(30) 



dt 



The next step is to determine an appropriate (approximate) representation of the ratio ip{^,^)/(f>{^,^). 
This consists in analyzing in details the eigenvalue problem for the adjoint operator of the linearization 

at{U{x,^),v,). 

First eigenfunction of the adjoint operator. As a first step, let us derive the specific expression of 
the hnearized problem. Given W{x,^) := (?7(x, ^), v*), let us look for solution to ([TT] ) in the form 
W + (u, v). The perturbation (u, v) satisfies 
du _ dv dU d^ 
dx dt ' 



dt 
It 



dx 



where 



d ( du dv 

— p(x) + a\{x) u + a2{x) v + b\{x)— — I- b2{x)— — I- h.o.t. 
dx^ dx 

p{x) 
a\(x) 
a2{x) 
bi{x) 
b2{x) 



= -F{U{x),v,) - E{G{U{x),v,)y 

- -daF{U{x),v,) - sduuG{U{x),v,)U'{x) 

- -dvF{U{x),v,) - EduvG{U{x),V,)U'{x) 

- -EduG{U{x),v*) 
= EV{U{X)) U~\x) 



and h.o.t. collects higher order terms. Then, the linear operator is 



-1 

dx \\a\{x) a2{x) 











b]_{x) b2{x) 



Setting w = {u, v) and taking the scalar product against the line vector a> - {(p, i/'), we get 



where 



- -r3'\w, + {£*m) ■ w 
dx ^ 



-1 

i\ a2 I \v 



b2 



(<!>' <A') 



0\/m 

^2 V 
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/O -ai(x)\/<^'\ d j(0 biix) 



and 

The Unearization of the boundary conditions ([TSl l for the operator gives 

(31) (f/v - v*m)(±0 = 0, v(-^) = 0. 

Thus, there holds 

J[w,co]l^_^ - (A(-a^i"' +^2V'||^_^ = 0, 

wj^.^^^ - -(P{+{M+0 + il>i+£){aiu + a2V + b^u' + ^2v'||^^^^ = 
and the requirement J'\^j = is satisfied if 

(32) (Pi+i) = 0, (A(±0 - 0, 

that are the boundary conditions for £.*. 

Next, we consider the eigenvalue problem for the adjoint operator with the aim of determining 
an approximate expression for the ratio i/^(^)/0(^), where (0, ip) is the eigenvector relative to the first 
eigenvalue A. The eigenvalue problem reads as 

( - ai{x)ilj' + {bx{x)i]j')' = A(P, 

1 0' - a2{x)tfj' + {b2{x)ilj')' = Alt/. 

with boundary condition (l32l ). Setting := tp' and w = {<p,tfj,G), after some standard algebraic 
manipulation, the above system can be rewritten as 

' AUIv, A an 
1 
-E-^AIduG s-^duFlduG^ 



dw 

(33) — ^Aw where A := 

dx 



with duF = duFiU, V,), duG = duGiU, V,) and 



v^,\U'^ j U dx 

In the limiting regime e — > 0^, the profile U tends to a step profile joining atx = ^ the values u±. Thus, 
we consider the matrix- valued A as piecewise constant and a representation for the eigenfunction can 
be obtained by an appropriate matching at x = ^ of solutions to (1331 ) with U considered as a constant. 

In order to determine the spectral decomposition of A and its behavior in the regime e ^ 0, we 
analyze the characteristic roots fi, solutions to 

. „ , (duF AU\ o 2AU A^ 

\<9„G V* / v{U) duG 

Plugging expansion for fj. in the form fi = yUie~' + /J-q + o(1) as e and skipping the calculations 
for shortness, we determine asymptotic formulas for the three roots jUo,iu± 

\ d F 

(34) ^ = = - ^— + o{£~^), ju = iu± = C±A + oil). 

£ OuG 
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with C+ are such that 

duF{U,v,)C^ + 2v, UC-\=Q, 
that gives, thanks to the specific form of the function F, 

2 



(35) 



Corresponding expression for the right/left eigenvectors r/^of A relative to the eigenvalues yUo,± can 
be determined starting from the relations 



- rs ^ 0, 
ArY + (sduGfi - duF)rj, = 0. 



AU \ I A 
yu n - - TT^rj ^ 0. 



£ duG 



Ari-idr2=0, 

where r = {r\,r2, ra) and € = {£1,^2, £3)- Taking advantage of the expansions (l34l ). we infer for /jq the 
expression 

ro = (0, 0, 1) + o(l), 4 - i-A/duF, 0, 1) + o(l) 

and for //+, using the notation given in (l35l) . 

1 



:(C±5„F, 1,C±^) + o(l), 4 - (C±, 1,0) + o(l) 



" 1 + C25„F 

Thus, disregarding the o(l)-terms in the expression for the right and left eigenvectors, we may intro- 
duce the projection matrices 

^ 


-A duF) 

^ jClduF C^duF 0\ 

r+ ~ 



1 



C± 
I ClA 



1 

C+/1 



" 1 + ClduF 

Finally, we get an approximated expression for the exponential matrix of A 

Considering the solution in the sub-interval {-£, ^) and taking into account the boundary condition 
il^i-i) = 0, we deduce 

Hence, we deduce 



m 
m 



duF 



Cl 



1 + ClduF 



cl 



1 + ClduF 



c 



m 1 + ClduF 1 + ClduF 

where 5„F is calculated at (m-, v*). For A ^ 0, there holds 

<A(^) C. . C+ , (C- + C+)(1+C_C+5„F) 



-I- 



0(0) 1 + ClduF 1 + C25„F 



+ o(l) 



(1 + Ci5„F)(l + ClduF) 



+ o(l)-o(l) 
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having used the relation C C+ = -1/(9„F. Thus, we infer from the expansion for fi±, the relation 

m i ci , ci m i^+i) , , 

w < ^ :i f- (c + + o{a) = A + o(A) 

Hence, we end up with the asymptotic expression 

¥S ^ ^ + i , 

<p{0 duF{u^,v^) 
Plugging into ( [30l ). we deduce the new form 

^3^^ d§_ ^ {^ + i)AmK^{^)-K-{^) 

dt duF{u-,Vt) u+ - u- 

where Ai{^) denotes the first eigenvalue of the operator 

Note that, since di,F{u-,Vt) < and - /<-_ is an increasing function, the equilibrium point ^* is 

stable if and only if Ai{^) < 0, as it should be. 

Further simplification of tlie slow motion equation. As a last step, we derive an approximate 
expression for the functions k± showing that the motion is indeed exponentially slow. Such functions 
are implicitly defined by the relations 



f 



duG{s,v,) , ^ + £ 
■ as - -- 



III, K+-F{S,V^:) £ 

In the limit e — > 0^, the righthand side blows up, hence k± is such that 

lim K+{^) = F{u+,Vt). 

£-->0+ 

To get a preciser description of k± for small e, we approximate functions F and duG by 

F(s,Vh.) w F(m+,v*) + duF(u+,Vt)(s - u+), 
duG{s,Vt) X duG{u+,v^) 
Setting h+{^) :- k+(^) - v*), we obtain 



•Jiu 



duG(u+,Vt) , ^ + t 

■ CIS = 



h±{^) - duF{u+, v^){s - u+) 

Integrating, it follows 

duG{u+,Vt) ( duF{u+,v^){u^ - u+)\ ^ + 1 
In 1 - 



from which we infer 



, ... duF (m, - U-) duF+{u+ - u,.) 



where = v*) and = duG{u+,Vt). Finally, since f (m_, v*) - v*), we obtain 

K+{^)-K-{^) X duF^{u+ - M*)exp-^-— — \ + duF (m* - u-)&x^\——s 

I duG+ e I \duG- e 
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as £ — > 0. Inserting into (l36l ). we obtain the equation 

(37) — - ^^V\-^-p^ + ^^P (^ + 0/li(^) 

which is th ultimate version of the equation for the dynamics along the approximate manifolds. 

4. Conclusions 

Equation (|37] ) describes the motion's equation along the approximate manifold {iy(-,^)) and it 
should be considered as an equation for the movement of the transition layer from v,.) to v,.) 
in the bounded interval {-i, £) due to the boundary conditions. To our knowledge, this is the first 
attempt of determining a relation that should be capable to quantify the boundary effects in the case 
of isentropic Navier-Stokes equation for compressible fluids. Starting from it, a number of possible 
issues amenable to future investigations arise. We propose a list of eventual topics, dictated by our 
personal taste. 

1. The determination of the size and, mainly, the sign of the first eigenvalue Ai of the linearized 
operator has not yet been explored. In particular, its sign is fundamental to determine stability 
of the single steady state. Comparing with the case of scalar Burgers equation, it is expected that 
asymptotic stability holds and thus that Ai is negative. 

2. The manifold builded to determine the form of the motion's equation is approximated. The 
residual is in any case exponentially small and thus the dynamics remain confinated in a small neigh- 
borhood around such manifold. In analogy with the results proved for the Burgers equation, it would 
be interesting to analyze if there exists an exact invariant manifold, lying close to the manifold of the 
approximate solutions, and determine a corresponding equation for the dynamics along it. 

3. The equation for the location f of the transition layer is almost explicit. The unique term 
remaining unknown is the first eigenvalue Ai that could be determined numerically. We wonder if it 
would be possible to verify numerically the reliability of the equation, bypassing the smallness of the 
right hand side and, thus, the very long timescale of the shifting phenomenon. 

4. As shown at the beginning of Section 2, the approach is formal but flexible and, in principle, ca- 
pable to be applied to different systems fitting into the general class of hyperbolic-parabolic systems. 
Among others, the extension of the approach to the case of non-isentropic Navier-Stokes would be 
particularly intruiguing. 

5. Following results available in the literature (see Introduction), the case of the half-line could 
also be taken into account. In particular, to quantify by means of a motion's equation the effect of the 
presence of a single boundary seems to be a natural variant of the case considered here 
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